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A MAXIMUM  LIKELIHOOD  KST 1 MATOR  FOR  AN  EXPONENTIAL 


PARAMETER  FROM  A LIFE  TEST  WITH  BOTH 
TYPE  I AND  TYPE  II  CENSORING 

Abstract 

A hybrid  life  test  on  items  assumed  to  have  an  exponential 
lifetime  combines  type  1 and  type  II  censoring.  In  type  I 
censoring,  n items  are  placed  on  test  and  observed  for  a fixed 
time  t*,  while  in  type  II  censoring  the  test  terminates  with 

the  r^*1  failure,  where  r is  a preassigned  integer.  If  tf 

th  ° ° 

is  the  time  of  the  rQ  failure,  a hybrid  life  test  terminates 

at  min(t  ,t*).  In  some  situations  it  may  be  of  interest  to 
ro 

estimate  the  average  lifetime  G , following  the  test  deci- 
sion. In  this  report,  we  find  the  maximum  likelihood  estima- 
tor, 0 , when  the  sample  is  subject  to  hybrid  censoring. 

An  expression  for  E(Q)  is  derived.  Because  this  expression 
is  complex,  computer  simulations  are  used  to  examine  the  bias 
of  0 . 


A MAXIMUM  LIKELIHOOD  ESTIMATOR  FOR  AN  EXPONENTIAL 
PARAMETER  FROM  A LIFE  TEST  WITH  BOTH 
TYPE  1 AND  TYPE  II  CENSORING 


I . Introduction 


Collecting  data  for  a life  test  is  often  complicated 
by  some  type  of  censoring  on  the  observed  lifetimes.  Con- 
sequently, testing  schemes  involving  various  combinations  of 
time  and  failure  censoring  have  been  developed.  In  type  I 
censoring,  n items  are  placed  on  test  and  observed  for  failures 

for  a fixed  time  period  t*.  In  type  II  censoring,  n items 

^ h 

are  placed  on  test  and  observed  until  the  r failure  occurs, 

o 

where  r ^ n.  The  value  of  r is  chosen  prior  to  the  test, 
o o 1 

Epstein  (1954)  proposed  a testing  scheme  which  combines  type  I 

and  type  II  censoring.  In  his  scheme  the  test  terminates 

at  min(t  ,t*),  where  t is  the  time  of  the  r ^ failure, 
r r o 

o o 

We  shall  refer  to  this  testing  scheme  as  the  hybrid  test. 

In  this  report  we  shall  assume  that  the  lifetimes  under 
consideration  have  an  exponential  probability  distribution. 
Thus,  if  : represents  the  lifetime  of  an  item,  the  density 
function  of  \ is 


f 


f (t; 0) 


l - t/e  . * 

— e i f t 

0 


(l.D 


elsewhere 


where  0 > 0.  In  this  form  9 is  the  average  lifetime  of 


an  item  in  the  population. 
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Epstein  (1954)  developed  the  hybrid  scheme  to  test  Ho:  0 = Gq 

ciqainst  II  ^ i»  - G ^ , where'  0 < 9 . If  min  (t  , t*)  t*  , 

o 

H is  accepted,  while  if  min(t  ,t*)  = t ,11  is  rejected, 
o r r o 

° o 

figure  1.1.  illustrates  this  hybrid  testing  scheme. 


fail ures 


Figure  1.1.  Possible  sample  paths  in  a hybrid 
scheme . 

It  is  very  likely  that  a true  value  of  0 less  than 
may  lead  to  rejection  of  Hq,  or  a true  value  of  0 greater 
than  0q  may  lead  to  the  acceptance  of  Hq.  Consequently , 
it  may  be  of  interest  to  compute  an  estimate  of  9 following 
the  life  test  decision  using  the  test  data.  The  estimate 
might  be  a point  estimate  or  an  interval  estimate  of  9 . 

In  this  report  we  shall  consider  the  problem  of  finding 
a point  estimator,  the  maximum  likelihood  estimator  of  9 , 
when  the  data  is  collected  for  the  hybrid  life  test.  We 
shall  assume  that  the  life  test  is  conducted  without  replace- 
ment, i.e.  failed  items  are  not  replaced  with  new  items. 
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Point  estimation  for  0 under  type  J and  type  II  cen- 
soring is  thorouyhly  discussed  in  the  literature  of  reliability 
and  life  testimj.  The  maximum  likelihood  estimator  of 
under  type  II  censoring  was  found  by  Epstein  and  Sobel  (1953) . 

If  t,,t„,...,t  represent,  the  ordered  failure  times  of  the 
12  n 

n items  on  test,  the  maximum  likelihood  estimator  of  0 is 
given  by 


1 ° 

0 = - l I!  t. 
r . , i 

o 1 = 1 


(n  - r )t  ] = T /r 
o r r o 

o o 


(1.2) 


where  T 

r 

o 


r 

o 


l t . 

l 


i=l 


+ (n 


r ) t represents  the  total 
o'  r 1 
o 


accumulated  test  time,  or  total  time  on  test,  at  the  time 

of  the  r ^ failure,  t . The  maximum  likelihood  estimator 
° 

of  0 under  type  I censoring,  and  its  properties,  are  dis- 
cussed by  Bartholomew  (1957),  Mendenhall  and  Lehman  (1960) 

and  Bartholomew  (1963),  including  a more  general  situation 

* 

where  each  item  on  test  has  its  own  truncation  time  t^  , 
i=l,2,...,n.  Lnder  type  I censoring  the  maximum  likeli- 
hood estimator  is  given  by 

1 k 

6 = i[  r t.  + (n  - k)t*]  (1.3) 

K i=l 


where  k is  the  number  of  failures  observed  by  time  t*. 
These  results  and  properties  of  the  estimators  are  sum- 
marized in  Mann,  Schafer,  and  Si ngpurwalla  (1974). 
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We  may  conjecture  that  (1.2)  and  (1.3)  could  be  combined 
to  form  a maximum  likelihood  estimator  under  the  hybrid 
censoring.  This  conjecture  will  be  shown  to  be  true  in 
this  report.  For  each  item  on  test  under  this  hybrid  scheme, 
a mixed  distribution  is  defined  having  positive  probability 
at  the  point  t* . Order  statistics  from  this  distribution, 
which  admit  the  possibility  of  ties,  are  defined  and  the  like- 
lihood of  the  first  r order  statistics  is  found.  An  ex- 

o 

a. 

pression  for  E(0)  is  also  derived  in  section  II. 2.  A simple, 

closed  expression  for  the  bias  is  found  when  r = 1 . A 

o 

A 

closed  but  complex  expression  for  E(6)  is  given  for 

2 < r ^ 5.  The  bias  of  0 is  further  examined  by  simu- 
o 

lation,  using  some  of  the  actual  hybrid  test  schemes  of 
Epstein  (1954) . 


II.  A MAXIMUM  LIKELIHOOD  ESTIMATOR 


FOR  THE  HYBRID  Tf 


II. 1.  Derivation  of  the  maximum  Likelihood  Estimator 


under  a hybrid  test  which  admits  tics  with  positive 


statistics  is  then  found,  from  which  the  maximum  likelihood 


Now  define  a measure  y as  follows:  If  L is  Lebesque  measure 
and  n is  the  counting  measure,  then  for  any  Lebesque 
measurable  set  A, 


Then,  if  we  define 


6 


f(t)  is  a density  in  the  sense  that  it  is  the  Radon-N ikodyn 

derivative  of  P with  respect  to  p.  If  i _ , i _ , ....  t are 

1 2 n 

he  n unordered  failure  times  of  the  items  in  the  hybrid 
test,  we  may  consider  . ..,  to  be  a random  sample  from 

a distribution  with  density  function  (2.1).  Then  the 
likelihood  of  6 is  given  by 


L ( 0 ) = n f . (t . | G) 
i = l 1 1 


Let  0 < t,  < t_  •••  < t < t*  define  the 

corresponding  order  statistics.  We  have  the  possibility 

of  ties  here  since  an  item  may  assume  the  value  t*  (i.e., 

-t* /0 

censored  lifetimes)  with  probability  e . Consequently, 

the  problem  of  these  possible  ties  must  be  considered  when 

finding  the  joint  density  of  the  order  statistics.  In  the 

case  of  ties,  the  order  statistics  are  defined  uniquely, 

though  somewhat  artificially,  in  the  following  manner.  If 

and  t are  tied  at  t*  for  the  kth  ordered  failure  time 

and  i < j,  define  t.  = t.  and  t,  , , = t . . In  most 
J k l k+1  3 

considerations  of  order  statistics,  the  density  is  assumed 
to  be  continuous,  and  points  in  n-space  where  two  or  more 
ties  occur  are  ignored  since  they  have  zero  probability. 

Now  these  points  must  be  considered  separately.  To  that 
end,  we  partition  n-space  into  sets  AQ , A^  , ...,  A^,  where 
we  define  A^  = {t^,  •••»  ’ exactly  i points  equal  t*}. 
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Thus,  Ag  is  the  set  where  all  coordinates  differ  from  t*  , 
while  A-j,  for  example,  is  the  set  of  points  where  three 
coordinates  equal  t*.  The  joint  density  of  the  order 
statistics  will  differ  on  each  A^.  We  further  partition 
each  A^  to  obtain  the  required  one-to-one  transformation 
onto  the  set 


l(tl'  L2 V ; ° < fcl  - fc2  - - Si  - tM 


For  example,  when  n = 3,  A^  is  partitioned  into 


a21  = {(V  V T3);  0 " 'l  < T2  = T3  = tM 

* = Tl'  fc2  = V fc3  = T 3 ' 

a22  = {(V  T 2 ' T3);  ° " T2  " T1  = T3  = tM 
^ tl  ~ X2'  t2  = Tl'  t3  = T3' 


and 


a23  ‘ {<V  T2 ' T3>  ! 0 ‘ t3  ■=  T1  * t2  * tM 

^ ll  1 3 ' l2  ' V C 3 = 1 2 ’ 


Consequently,  on  A^  with  n = 3,  the  joint  density  of  the 
order  statistics  is 


cj(t1,  t2,  t3)  = 3f  (t1)  f (t2)  f (t3)  . 
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In  general,  there  will  be  (^) (n-i) ! sets  in  the 

partition  of  to  accomplish  the  one-to-one  transformation. 

(^)  accounts  for  the  number  of  ways  to  select  the  i ties, 

while  the  (n-i)  untied  coordinates  can  be  permuted  (n-i) ! 

ways.  The  joint  density  for  the  order  statistics  over  set 

A.  can  be  written  as 
1 


g(t1/ 


V 


(")  (n-i)  1 


n 

n f(t.) , i 
j=i  3 


Of  Xf  • • • f n • 


If  we  define  C . = ( *?)  (n-i)  i = ni/i!  and  let  q = e 
n f jl  x 

then  the  joint  density  of  the  order  statistics,  in 
generality,  is  written  as 


-t*/Q 


full 


i 


n 


c n n f(t.)  if  o<t, <t  <• • • <t  < t* 

n ,0  l 1 2 n 


n-1 

C , q n f(t.)  if  0<t. <t  <• • • <t  = t* 
n'x  i=X  1 1 2 n 


n-2 


' _ q II  f(t.)  if  0 <t.  <t_<  • • • <t_  , = t 

n,2  ^ i 12  n-1  n 


= t* 


C q 
n,n^ 


n 


if  0 <t.  = t_  = •••  = t = t*. 
x 2 n 


v - 
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It  is  now  possible  to  find  the  joint  marginal 
density  of  the  first  out  of  n order  statistics  from 
density  ( 2.1  ).  Recall  that  r^  is  the  truncation  value  on 

the  number  of  failures  for  the  hybrid  test.  This  joint 
marginal  density  can  bo  found  from  the  theorem  which  follows. 

Theorem  2.1.  If  f(t)  represents  the  density 
function  (2.1  ) and 


t 

F ( t)  = /f(t)du  , 
0 


The  last  equality  follows  from  the  observation  that 


i l 


0 < 4 * 


t , = t * t*  , 
n- 1 n 


the  integration  will  be  over  the  single  point  t*.  Thus, 


g(  t , . . . , t )=C  / fl  f(t.)q  dii 

i ni  n ' ^ { t*  } i = 1 1 


= cn  2 q n f(ti) 
n'^  i=l  1 


if  0 < t.  < t _<•••<  t , = t * 
12  n- 1 


Likewise , 


g(tl tn-l)  = Cn,3q  i!l1f(ti)  if  0<tl<t2<,,,<tn-2=tn-l  = t* 


g(tl' * * * 'tn-l)  = Cn,nq 


if  0<t1=t2=---=tn_1=t*, 


which  is  the  desired  result  for  r = 1.  Now  assume  the 
theorem  is  true  for  r = k.  It  suffices  to  show  it  also 
holds  for  r = k + 1.  Again  in  the  region 


0 < t.  < • • • < t , < t* , 

1 n-k  — 


9(tl'  * " ' tn-k-l) 


n-k- 1 t*- 

= c n f ( t ) / f 

' i-l  t . = t . , 

n-k  n-k- 1 


f(t  . ) [ 1- F ( t . ) Pdt  . 
n-k  n-k  n-k 


Then  in  the  region  whore  0 < t,  < •••  < t , , = t , t* 

1 n-k-1  n-k  ' 


g(tl*  V-k-l*  ~ i(*]Cn,k+2 


n-k-2 


i = 1 


f . . . k-t  2 , 

f(t.)q  du 


k + 2 

Cn,k+2  g if1  f<V 


if  0 < t,  <•••  <t  . . = t * . 

1 n-k-1 


Likewise , 


9(tl Vk-l*  ‘ Cn,k+  3 q 


k + 3 


n-k-  3 
ii  fit.) 
i = l 1 


if  0 < t < •••  < t . = t . = t* 

1 n-k-2  n-k-1 


g<tx. 


t . , ) = C q 

n-k-1  n,n 


if  0 < tj  = t2  = 


'"n-k- 1 


t*  . 


This  establishes  the  result  for  r = k + 1 and  the  proof  is 
complete  by  induction. 

From  Theorem  2.1  , the  joint  marginal  density  of 

the  first  r^  order  statistics  is  seen  to  be 


15 


whore  k = the  number  of  failures  observed  in  (0,t*).  Note 
that  0 < k < r and  k = r if  and  only  if  t < t*. 

- 0 0 ro 

Consequently  we  may  write  the  likelihood  function  of  the 
first  r order  statistics  as 

k a(n-rQ) 

UO)  « I II  f (t.  ) ]q{1_a)  (n'k)  [1-K(t  )]  (22 

i = l 1 *o 

where 


a 


1 

< 

0 

i 


if  t < t* 
r0 

if  t > t*  . 
0 


In  terms  of  density  function  ( 2.2  ),  the  likelihood  becomes 


1(0) 


k 


-~(l-a)  (n-k) 
e 


( n-r. 


) a 


k 

l l t . + ( 1‘ a)  (n-k)  t*+a  (n-r. ) t ) 
i=l  1 0 ro 


Then , 


' i VS*  + 


(l-a)(n-k)t*  + a(n-rQ)tr  1 


lnt (G)  = C - klnO 


and 


]6 


BlnL(O)  _ -k 
39  6 


k 

I l t + ( 1-a) (n-k) t*+a (n-r„ ) t 
1 Or 


] 


Setting 


31nL ( Q)  _ 

30  " u' 


we  find  that  the  maximum  likelihood  estimator  for  6 is 

k 


6 - I t + ( 1-a)  { n-k)  t*  + a(n-r„)t 

0 r0 


. ( 2.3  ) 


This  is  merely  the  combination  of  the  maximum  likelihood 
estimators  from  the  type  I and  type  II  censoring  which  we 
expected . 

11,2  The  First  Moment  of  the  Maximum  Likelihood 
Es timator  ~ — ~ 

The  maximum  likelihood  estimator  derived  in  the 
previous  section  can  be  used  to  estimate  0 following  the 
hybrid  life  test.  In  order  to  examine  the  bias  of  this 
estimator,  an  expression  is  derived  in  this  section  for  the 
first  moment  of  Q.  In  section  II. 1 we  obtained 

a 1 k 

6 = k * 1 + (1-a) (n-k) t*  + a ( n- r_ ) t }, 

i=l  0 rn 
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f • • • / 


whero  ti  = ith  ordered  failure  time,  i - 
t*  = truncation  value  for  time, 
r = truncation  value  for  failures,  and 
k0  = number  of  failures  at  time  of  decision. 

The  distribution  of  k is  given  by 


n. 


P(k  = j)  = <1 


(")pj(l-p)n'j 


if  j = 0,  1,  •••»  r0'1 


. Z (")PXU-P)n'  if  3 " r0' 

{ 1=ro 

where  p = P(an  item  fails  before  t*)  - 1 G 

(23),  "a"  is  a random  indicator  with  a distribution 

\ - D , v _ r \ = pit  <t*)  , and  P(a=0)  = 

given  by  P(a-l)  - P(k  rQ)  rQ 

1-P  (a=l)  . 

The  estimator  9 will  be  infinte  when  k = 0,  so 
we  will  condition  on  5 being  finite  and  find  E(G|S<»>  = 
E(0 |k>0) . Since  e is  infinite  with  probability  f(k-O). 
the  distribution  of  5 conditional  on  6 < » is  tho 

unconditional  distribution  of  5 divided  by  the  constant 
POl>0).  Hence,  E < 5 | 6 <-)  will  be  merely  the  unconditional 
expectation  of  § (over  finite  values  of  6)  divided  by  the 
constant  P(k>0).  So,  E(S|G<<») 


P(k>0) 


E E t. 


a (n-r  ) t 

, t.  (1-a)  (n-k)  t*  . „ u r0 

+ L k + E k 


( 2.4  ) 


Each  of  the  expectations  in  ( 2.4  ) will  be 

evaluated  in  turn.  We  first  write 


k 

E E t 


0 j 

= l E l t k=j  P (k= j) 

L*-H  J 


Given  that  k = j,  for  j = 1,  2,  ...»  rQ-l,  there  are 
exactly  (^)  equally  likely  ways  of  selecting  the  set  of  j 
items  out  of  n which  will  fail  first.  If  we  condition 
again  on  each  of  the  (?)  equally  likely  combinations  it  is 
seen  that 


E E t |k=j  = E | E T k=j  = E E T.  , 

i-i-rl  U-i  41  i-i  -4 


where  is  the  i unordered  failure  time  whose  distribu- 
tion is  truncated  at  t*.  The  truncated  distribution 
follows  from  the  fact  that  we  are  given  the  number  of 
failures  observed  before  time  t*  is  k = j.  Thus  the 
distribution  of  is  given  by 


Direct  integration  using  ( 2.5  ) gives  the  result  that 


Now,  k = r implies  that  at  least  r.  items  failed 


by  time  t*  (even  though  observation  ceases  at  t ) . We 


must  treat  this  situation  differently 


To  evalute  this  expectation  we  need  the  joint  density  of 
the  first  rQ  order  statistics  taken  from  the  exponential 
population,  conditioned  on  t < t*.  Let  x.  represent 


the  observed  value  of  the  random  time  t 


Then 


< min ( xr. , t*) ] 
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t* 


rQP(tr  <t*) 


t* 


t* 


( l x.) 
i=l 


n ! 


Xj=0  x2=xx 


X =x  . 

ro  ro'1 


(n-rQ)  10 


e 


1 

0 


r 


0 


l F. 

i=l 


x^  + (n-r 


0 


( 2.6  ) 


Evaluating  and  simplifying  this  r^-fold  integral  becomes 
very  tedious  as  increases.  We  have  evaluated  ( 2.6  ) 
for  values  of  r^  through  five.  We  give  here  the  results 
for  r^  = 1 and  2,  while  the  results  for  r^  = 3,  4,  and  5 are 
given  in  Appendix  1.  In  these  expressions  we  define 


C = nj . ( 2.7  ) 

r (n-r)  ! rP  ( t <t*) 

The  integral  ( 2.6  ) reduces  to 


nt* 

f 1 

e 

0+nt* 

9 

2 

2 

e 

n 

l n J 



when  r^  = 1 , and 


0 ( 3n- 2 ) 

r 

2nt*+  (n  + 2)  0 

nt* 

e 6 - 

] 

n0+  (n-1) t* 

(n-1) t* 

9 

2 2 
n (n-1)  * 

2 

l n 

C 

. (n-1)2 

G 

We  now  have,  for  the  first  expectation  in  expression  (2.4  ) 


k 

Z t. 

ro_1 

= z -(-°~2t.*)  P(k=i)  + E 

k 

IT1  irj 

i=l  p 

b irl  °J 

P(k=rQ) 


P (0  <k<rQ)  (0-qt*) 


+ E 


E t.  |, 


P (k=rQ) 


( 2.8  ) 


where 


k 

Z t. 

i=l  j~l  k=rQ 


is  given  by  (2.6  ) . 

Next,  we  evaluate 


{ 


(1-a) (n-k)  t* 
k 


3- 


If  we  condition  on  the  value  of  a,  this  expectation  becomes 


(r 


E r-*-1  a=0 


P (a=0)  = 


E[¥|k-0>*'V  - t-PO'-'r,,) 


'o'1 


nt*  Z (f)P(k=i)  " t*P (k<r„)  = 
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i 


VJ 

t*  z -(-n~I)  p(k=i; 
i=i 


To  establish  this,  we  made  use  of  the  fact  that  a = 0 if 
and  only  if  k < r . As  before. 


t* 


t* 


P(k---i)  = (")  ( 1-e  0 ) 1 (e  0 )n  1 


Finally,  we  find 


at 


<n-r0)  D 


Conditioning  again  on  the  value  of  a,  we  get 


rat  n 

t 

t 

ro 

= F 

r°  | a--l 

P (a=l ) = E 

ro 

k 

Ci 

k * I 

r 

L J 

0 

“| 

t_  <t*  i p ( t <t*)  , 
0 


since  a - 1 -*•  t < t*  ■+  k = r . Here  we  require  the 

0 u 

density  of  the  order  statistic  from  the  exponential 

distribution  which  is  truncated  at  t*. 


f (x  ; |t  <t*)  = 
r0  0 0 


r 

X 

ro 

"O  1 (n-rn  +l)x 

u r0 

n ! 

1-e  0 

e 0 

P(tr^<  t*) 0(n-ro)  ! (rQ-l)  ! 


if  x < t* 


0 elsewhere. 
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Combining  the  three  parts  from  ( 2.4  ) we  have  a final 

A ^ 

expx'ession  for  the  first  moment  of  0,  conditional  on  G 
being  finite. 


A „ P (0<k<rn) (Q-qt*)  P (k~rn ) 

E ( 0 | G<°°)  + E 


P{k>0)p 


rQP(k>0) 


E t.  t < t* 
i=l  1 r0 


rn_1 

+ fc- 

0 

Z P 
i=l 

P (k  >0 ) 

n-1 

(n-rQ)  n9 

V1 

rQ  P (k>0 ) 


r -1  . |r0-1'' 

? 1-1)  j 1.3 


j=0  (n-r0+j+l) 


_ (n-rQ+ j+1) t* 


1-e 


(n-rQ+ j+1) t* 


+ 1 


( 2.10) 


In  the  simplest  state,  when  r^  = 1,  E ( 0 | k >0 ) 
reduces  to 


where  P(k>0) 
So , 


__1 

P(k>0) 


P(t  , 


te-(e+nt*)e“nt*/0]  , 

i -nt*/G  . 

< t* ) = 1-e  = p* , 


say. 
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When  r = lf  0 has  a negative  bias  equal  to  nt*  (^*-  -1) . 
Similar  simplifications  when  > 1 were  not  found. 

If  this  maximum  likelihood  estimate  of  0 is 
computed  with  the  data  from  a hybrid  life  test,  the  values 
of  n,  Tq , and  t*  are  supplied  by  the  test  and  cannot  be 
adjusted  for  the  purpose  of  reducing  the  bias.  We  have 
simulated  some  hybrid  life  tests  on  the  computer  for  the 
purpose  of  observing  some  bias  values  which  the  maximum 
likelihood  estimator  inherits  from  the  test.  I f we  test 
Hq  : 0 = 9g  against  : 0 = 0^,  where  0^  < 0q,  Epstein 
(1954)  has  given  the  theory  necessary  to  find  the  values 
of  n,  Tq  and  t*  which  correspond  to  given  error  probabili- 
ties a and  B.  Simulations,  of  1000  trials  each,  were  run 
on  27  different  hybrid  test  schemes.  In  each  simulation 
E ( 0)  was  estimated  and  the  bias  was  estimated  at  0 = 0^. 
Table  2.1  gives  the  simulation  results.  The  bias  is 
expressed  in  the  table  as  the  percentage  of  the  actual 
value  of  0.  Using  1000  trials,  the  estimated  standard 
deviation  for  this  percentage  was  less  than  3%  for  all 
simulations.  For  comparison,  the  table  also  gives  the 
exact  bias  in  each  case  for  the  maximum  likelihood 
estimate  from  a type  I testing  scheme  having  the  same 
values  of  n and  t*.  This  exact  bias  was  shown  by 
Bartholomew  (1957)  to  be 
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nt*E  ( j~)  - 3—  - t* 

K p 

whore  k = the  number  of  failures  in  (0,t*).  It  is  assumed 

that  k > 0 and  E(l/k)  is  conditional  on  k > 0 . As  before, 
-t*/0 

3 “ e an<3  P - l~q.  It  is  clear  from  Table  2.1  that 

the  bias  may  be  quite  large  especially  for  smaller  sample 


si zes . 


APPENDIX  1 


Listed  here  are  the  values  of  the  integral  ( 2.6  ) 


when  r 


equation  ( 2.7  ) 


2n0  + 30+ 3nt* 


20 ( 3n- 7n+  3) 
n2 (n-1) 2 (n-2) 


2n0-30+ (n-2) t* 
2 (n-2)  2 


2n0+  2 (n- 1) t 
(n-1)2 


3n0+4nt*+0 


20(5n-25n+35n-12) 


3n9-40+2 (n-2) t* 
2 (n-2) 2 


3n8+3 (n-1) t* 
2 (n-1) 2 


3n0- 80+ ( n- 3) t* 
6 (n- 3) 2 


(15n-130n+375n-404n+120) 
n2  (n-1)  2 (n-2)  2 (n-3)  2 (n-4)  2 


2n0+2 (n-1) t* 
3 (n-1) 2 


4n0+5nt*+50 
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